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PROPAGATION  OF  DISCONTINUITIES. ... 

IN  HETEROGENEOUS  ANISOTROPIC  PLATES 
By 

A.  S.  D,  Wang^  and  D.  L.  Tuckmantel2 
ABSTRACT 

Elastic  stress  waves  propagating  in  thin,  laminated  composite 
plates  are  analyzed  on  the  basis  of  a  lamination  theory.  The  theory 
is  based  on  the  Kirchhoff  assumptions,  but  it  includes  the  effects 
of  shear  deformation  and  rotary  inertia,  similar  to  Mindlin's  theory 
for  homogeneous  isotropic  plates.  The  individual  layers  comprising 
the  plate  are  assumed  to  possess  different  thicknesses  and  material 
properties.  In  particular,  each  layer  may  be  arbitrarily  anisotropic. 
Thus,  a  general  coupling  in  shear,  bending,  twisting  and  extensional 
effects  is  present  in  the  plate  constitutive  relations.  This  coupling 
results  in  simultaneously  coupled  stress  waves  propagating  in  the 
plane  of  the  plate.  Several  numerical  examples  involving  laminated 
fiber-reinforced  composite  plates  are  presented. 


■^Associate  Professor  of  Applied  Mechanics,  College  of  Engineering, 
Drexel  University,  Philadelphia,  Pennsylvania  19104. 
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I.  INTRODUCTION 


The  dynamics  of  heterogeneous  anisotropic  solids  has  been  an 
important  field  of  study  in  recent  years.  Increasing  popularity  in 
the  use  of  modem  composites  as  a  structural  material  has  necessitated 
intensive  research  into  their  material  characterization.  Their 
intrinsic  properties  are  being  investigated  from  both  a  phenomenological 
and  microscopic  point  of  view.  This  effort  has  resulted  in  the 
development  of  various  lamination  theories  describing  laminated 
anisotropic  plate  structures.  The  papers  by  Yang,  Norris  and 
Stavsky  [1],  and  Whitney  and  Pagano  [2]  are  among  the  most  notable 
works.  The  former  concerns  plates  made  of  layers  possessing  arbitrary 
anisotropy;  the  latter,  with  some  modifications*,  is  valid  for  plates 
made  of  monoclinic  layers,  i.e.,  each  layer  of  the  plate  possesses 
a  mid-plane  material  symmetry.  Both  of  these  theories  follow  the 
basic  approach  set  forth  by  Mindlin  [4]  for  homogeneous  isotropic 
plates,  and  include  transverse  shear  deformation  and  the  effect  of 
rotary  inertia.  Recently,  Chou  and  Carleone  [5]  improved  the 
assumptions  concerning  the  transverse  shear  deformation. 

Such  lamination  theories,  however,  describe  the  macro -character¬ 
istics  of  the  plate  rather  than  the  micro-characteristics  that  are 
effected  by  the  inhomogeneous  nature  of  the  composite.  Generally 
speaking,  the  macro-theory  is  accurate  for  plates  subjected  to  static 
loading,  but  is  deemed  inaccurate  when  applied  to  stress  wave  problems. 
This  is  true  in  so  far  as  the  wave  lengths  are  short  compared  with, 
say,  the  thickness  of  the  material  layers  of  the  plate.  But  for  low 
frequency  waves  with  length  longer  than,  say,  the  plate's  thickness, 
the  macro-theory  may  be  regarded  as  a  valid  basis  of  analysis. 

Based  on  this  assumption,  Moon  [6]  recently  investigated  stress 
waves  in  a  specially  laminated  fiber -reinforced  composite  plate, 
using  the  effective  modulus  theory  which  uncouples  the  transverse, 
bending  and  extensional  displacements.  A  wave  surface  approach  was 
used  to  describe  the  propagation  of  plane  acceleration  waves.  This 
method  yields  the  wave  velocity  surfaces  in  the  plane  of  the  plate. 


*A  comparison  of  the  two  theories  is  contained  in  Ref.  [3]. 
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The  present  paper  is  concerned  with  stress  waves  in  laminates 
that  are  composed  of  layers  possessing  arbitrary  anisotropy.  General 
shear,  bending,  twisting  and  extensional  coupling  is  present  in  the 
plate  constitutive  relations,  resulting  in  simultaneously  coupled 
wave  surfaces  in  the  plane  of  the  plate.  We  follow  general  laminated 
plate  theory  and  apply  a  control  volume  approach  for  the  analysis  [7]. 
Explicit  solutions  for  the  coupled  wave  surfaces  and  their  velocities 
are  obtained.  Several  numerical  problems  involving  laminated  fiber- 
reinforced  composite  plates  are  presented  and  their  unique  features 
discussed. 

II .  ANALYSIS 

Let  us  consider  a  thin  laminated  plate  of  thickness  h,  Figure  1. 
The  laminae  comprising  the  plate  are  assumed  to  be  individually 
homogeneous  and  anisotropic.  Thus  the  in-homogeneity  of  the  plate 
occurs  only  in  the  thickness  direction.  The  constitutive  relations 
for  any  one  of  the  laminae  are  given  by 

°i  -  Cij  ej  ^  "  1,2 >^>4>5,6.  Cl) 

where  the  C^.'s  are  the  elements  of  the  stiffness  matrix,  the  stresses, 

a ■  ,  are  defined  as  a,  =  a  ,  a0  =  a  ,  a,  =  a  ,  a.  =  a  ,  or  =  a  , 

x*  1  x’  2  y’  3  z’  4  yz*  5  zx’ 

a6  =  axy’  anc*  t^ie  strains ,  e . ,  are  defined  in  the  same  manner  as  the 
stresses . 

Following  the  theory  developed  by  Yang,  Norris  and  Stavsky  [1], 
we  assume  the  displacement  field, 


u  =  u°(x,y,t)  +  z  ipx(x,y,t) 

v  =  v°(x,y,t)  +  z  *  (x,y,t)  (2) 

w  =  w°(x,y,t) 

where  the  coordinate  system  (x,y,z)  is  shown  in  Figure  1  and  u,  v 
and  w  are  the  displacements  in  the  x,y  and  z  directions,  respectively, 
u°.  v°  and  w°  are  the  displacement  components  at  z  =  0,  in  the  x,y  and 
z  directions,  respectively  and  and  ijy  are  rotations  about  the  y 
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and  x  axis,  respectively. 

The  stress  and  moment  resultants  are  related  to  the  displacements 
by  (c.f.  Equation  14,  Ref.  [1]). 
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(Mx’  V  V  =  J  h/2  (ol*  °2  *  a6}  zdz 


(A...  B..,  D..) 
13  13  13 


f  +h/2  2 

C..  (1,  z,  zz)dz 

-h/2  J 


i,j  =  1,2,4, 5, 6  (5) 


and  the  notation  ip  ,  e.g.,  represents  partial  differentiation  of 

x>y  x 

with  respect  to  y. 

We  now  consider  a  wave  which  originates  at  an  arbitrary  point 
in  the  plate,  for  convenience  let  us  say  at  the  origin  of  the  (x,y,z) 
system,  and  propagates  in  the  x,y  plane.  At  any  given  instant,  the 
wave  surface  is  denoted  by  S,  as  shown  in  Figure  1.  Let  h  be  the 
normal  of  S  at  a  point,  A,  on  S,  and  let  s  be  the  tangent  of  S  at 
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the  same  point.  The  wave  surface  S  is  assumed  to  propagate  in  the 
direction  n  with  a  constant  speed  c. 

Let  u,*  v,  w,  u°,  v°,  w°,  i|»  and  be  the  displacements  (see 
Equations  (2)  ) referred  to  the  local  coordinates  (n,s,z).  The 
kinematic  (continuity)  conditions  across  S  at  any  given  point  on  S 
require  (see  [7]), 


[u°]  =  [v°]  =  [w°J  =  [*n]  =  [*s]  =  0 
tu’n’  v,n’  w,n’  *n,n»  ^s,irl  =  E 


^n,t’  ^s,t^ 


(6) 


and  (u,  v,  w,  vi>n,  4»s)  >s  =  0 

where  [  ]  represent  a  discontinuity  of  the  enclosed  quantity  across 

S. 

With  these  conditions,  the  plate  constitutive  relations, 
Equations  (3),  when  transformed  to  the  local  coordinates  (n,s,z) 
yield** 
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To  establish  the  dynamic  relations  across  the  wave  surface, 
let  us  define  a  control  volume  which  is  located  on  S  at  point  A, 

*Quantities  with  a  bar  on  top  are  referred  to  the  (n,s,z)  coordinate 
system. 

**For  the  transformation  of  A..,  B. .  and  D..  from  the  (x,y,z)  system 
to  the  (n,s,z)  system,  refer x^to  ^Ref.  ^  [8] . 
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as  shown  in  Figure  1.  Since  the  control  volume  moves  with  the  wave 
front,  an  observer  fixed  with  it  sees  a  normal  influx  of  mass  entering 
with  a  speed  Up  and  a  normal  efflux  of  mass  leaving  with  a  speed  Up 


where  the  subscripts  1  and  2  refer  to  the  properties  in  front  of  and 
behind  the  wave  front,  respectively.  Thus  the  steady-state  con¬ 
servation  of  mass  for  the  control  volume  yields, 


f  +h/2  rcn 

(p2U2  -  PlUx)dz  =  0  (9) 

•'  "h/2 

where  p  is  the  mass  density  of  the  material,  U  is  the  particle  velocity 
relative  to  the  wave  front  and  the  subscripts  1  and  2  refer  to 
properties  in  front  of  and  behind  the  wave,  respectively.  . 

It  is  noted  that  condition  (9)  is  satisfied  by  a  more  restrictive 
condition,  resulting  from  the  classical  thin  plate  assumptions, 
namely; 

p2U2  -  (10) 

The  force  and  moment  resultants  acting  upon  the  control  volume  must 
satisfy  the  equations  of  balanced  momenta.  These  are, 

tNnJ  =  |  +V2  CU2  -  Uj)dz  =  Pc2  [u,°]  *  Re2  ^n>nJ 
-h/2 

[Nns]  =  +h//2piui  (v2  -  vx)dz  =  Pc2  [v,°]  +  Rc2  [*s>nJ 
-h/2 

[Mn]  -  [  +V2  PlUx  CU2  -  Ux)zdz=  Rc2  [u,°]  +  Ic2  [^n>n]  (ID 

-h/2 

[M^]  =  +h^2plUl  ^2  "  ^ zdz  =  Rc2  +  1q2  ^s,nJ 

-h/2 

[Qn]  =  f  +V2  Piui  (w°  -  w°)dz  =  Pc2  Iw,£] 

-h/2 


9 


In  obtaining  the  above  relations,  we  have  retained  only  the  linear 
terms  of  the  displacements,  consistent  with  the  plate  theory.  In 
addition,  we  have  assumed  density  to  be  a  function  of  z  alone  in 
presenting  the  following  quantities 


r  +h/2  7  _  t  +h/2  2 

(P,  R,  I)  =  p,  (l,z,zZ)dz  =  P  (l,z,z  )dz  (12) 

J  -h/2  1  j  -h/2 


where  pQ  is  the  undisturbed  density  of  the  laminae. 

Using  the  relations  (7) ,  (8)  and  (10) ,  we  obtain  from 
Equations  (11)  a  system  of  five  linear  algebraic  equations  relating 
the  discontinuities  in  the  normal  derivatives  of  u°,  v°,  w°,  ij> 
and  \ps 
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In  order  for  a  non- trivial  solution  to  exist,  the  determinant 
of  the  coefficient  matrix  of  Equations  (13)  must-  vanish.  Thus,  five 
possible  wave  front  speeds,  c,  may  be  determined  for  any  given 
direction  n. 

Since  the  vanishing  determinant  represents  a  fifth  order  equation 

2 

in  c  ,  numerical,  rather  than  analytical,  techniques  must  be  used  to 
obtain  a  solution.  However,  it  is  of  interest  to  note  that,  if  the 
laminated  plate  has,  for  each  lamina,  a  monoclinic  symmetry  (i.e.  a 
plane  symmetry  with  respect  to  the  mid-plane  of  the  lamina)  the 
constants  =  A ^  =  B-^  =  =  0  for  all  directions  in  the  x,y- 

plane.  In  such  a  case,  [w,°]  is  uncoupled  from  the  system. 

Consequently, 

4  -  *ss/p  tl4) 
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where  c,-  represents  the  propagation  speed  of  the  discontinuity 
[w,°]  in  the  direction  ii. 

Furthermore ,  if,  in  addition,  the  plate  has  a  symmetry  with 
respect  to  the  x,y-plane,  the  bending -extensional  couplings  B^. 
become  zero  for  all  i  and  j .  Then,  Equations  (13)  separate  further 
yielding  two  quadratic  equations  whose  roots  are 


(15) 


(16) 

In  the  particular  case,  such  as  that  considered  by  Moon  [6] ,  the 
extensional  and  the  bending  rigidities  are  proportional,  i.e., 

^ij/®ij  =  =  constant»  i»j  =  1,2,6 

Equations  (14)  and  (15)  are  identical  causing  c^  ^  an<l  c3  4  to 
coincide  (see  Equations  (17)  and  (23),  Ref.  [6]). 
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III.  NUMERICAL  ILLUSTRATIONS 

For  the  numerical  illustrations,  plates  which  are  laminated  with 
unidirectional  fiber -reinforced  composite  layers  are  considered. 

The  material  properties  of  these  layers  are  described  by  the  following 
engineering  constants* : 

El  =  25  x  106  psi,  Ej.  =  106  psi,  GLT  =  0.5  x  106  psi 
VLT  =  ®*25,  VTT  =  0*35,  pQ  =  0.073  pci 


*These  values  are  typical  of  high  modulus  graphite -epoxy  composites, 
Such  material  layers  may  be  considered  as  being  square  symmetric. 
The  computation  for  the  plates’  rigidities  and  their  transformation 
to  local  coordinates  was  carried  out  following  the  outlines  in 
Refs.  [2]  and  [8]. 
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where  E  is  Young's  modulus,  G  is  the  shear  modulus,  v  is  Poisson's 
ratio,  and  the  subscripts  L  and  T  indicate  directions  parallel  and 
normal  to  the  fibers,  respectively.' 

We  have  considered  four  plates  each  of  which  is  made  of  four 
layers  having  different  lay-up  angles  and/or  lamination  sequences, 
namely:  a)  (079070790°) ,  b)  (079079070°), 
c)  (+307-307+307-30°)  and  d)  (+30°/-30°/-30°/+30°) ,  where  the 
angles  are  positive  when  measured  counterclockwise  from  the  x-axis 
to  the  fiber  direction. 

The  roots  of  the  vanishing  determinant  of  Equations  (13)  are 
determined  by  the  Newton  Raphson  technique.  In  general,  five  distinct 
roots  representing  five  possible  wave  speeds  exist  in  any  given 
direction.  However,  in  certain  preferred  directions,  the  material 
properties  are  such  that  repeated  roots  may  exist. 

Figures  2-5  show  the  wave  velocity  surfaces  in  the  first 
quadrant  of  the  x,y-plane,  for  (0°/90°/0°/90°) ,  (0°/90°/90°/0°) , 

(307-307307-30°)  and  (30°/-307-30730°)  laminates,  respectively. 

1/2 

The  velocity  is  non-dimensionalized  by  a  factor  of  (E^/pQ)  '  .  It 
is  pointed  out  that  the  slowest  velocity  surface,  which  is  associated 
with  the  transverse  displacement  w,  is  uncoupled  from  the  other  four 
surfaces,  since  the  material  is  monoclinic.  However,  all  the  other 
four  velocity  surfaces  may  be  severely  coupled.  On  each  of  the 
corresponding  wave  surfaces,  discontinuities  in  normal  forces,  shear 
forces  and  bending  and  twisting  moments  exist  simultaneously.  Their 
relative  magnitudes  may  be  determined  from  Equations  (13)  once  the 
wave  speeds,  c,  have  been  determined. 

Multiple  coupled  one -dimensional  stress  waves  in  a  heterogeneous 
plate  were  first  treated  by  Wang,  Chou  and  Rose  [9]  using  the  method 
of  characteristics.  Subsequent  experimental  investigations  in  the 
same  problem  has  not,  as  yet,  confirmed  the  multiple  wave  nature 
associated  with  laminated  plates.  Recent  experiments  conducted  at 
Drexel  University  using  a  low  frequency  ultrasonic  transducer . to 
produce  normal -to -the -plate  pulses  showed  only  two  distinct  -wave 
groups  traveling  in  the  plane  of  the  plate.  At  present,  a  definitive 
conclusion  cannot  be  drawn  due  to  the  limited  data  available. 
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Figure  2.  Wave  velocity  surfaces  for  the  C0°/90o/0°/90o)  plate. 
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Wave  velocity  surfaces  for  the  (0°/90o/90o/0°)  plate 
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Figure  5,  Wave  velocity  surfaces  for  the  (+30o/-30o/-30o/30°)  plate. 
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